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Abstract 

We investigate the spatially heterogeneous dynamics in the SPC/E model of 
water by using molecular dynamics simulations. We relate the average mass 
n* of mobile particle clusters to the diffusion constant and the configurational 
entropy. Hence, n* can be interpreted as the mass of the "cooperatively 
rearranging regions" that form the basis of the Adam-Gibbs theory of the dy- 
namics of supercooled liquids. Finally, we examine the time and temperature 
dependence of these transient clusters. 

More than thirty five years ago Adam and Gibbs (AG) proposed a theory to describe the 
dynamics of supercooled liquids [1,2]. In their approach they suggest that the system changes 
its configuration by the motion of independent "cooperatively rearranging regions" (CRR). 
Their main result is that the diffusion constant D is related to the temperature T and the 
configurational entropy of the system S con { by 

D ex exp(-A/TS coni ). (la) 

In the thermodynamic limit, they interpreted S^onf as ks log W c , where W c is the number of 
configurations accessible by the system and ks is the Boltzmann constant. More recently, 



W c has been interpreted as the number of basins in the potential energy surface (PES) 
accessible to the system in equilibrium, facilitating direct calculation of S con f by computer 
simulations [3]. The AG prediction has been tested and appears to be valid across a wide 
spectrum of liquids [4,5]. The AG theory also hypothesizes a relation between S^onf and the 
characteristic mass z of the CRR. However, the CRR are not precisely defined in the theory, 
and in the absence of a definition it has been challenging to test this aspect of AG theory. 

More recently, computer simulations on simple systems (such as Lennard- Jones mixtures) 
have shown that particles of high mobility tend to form clusters, and the concept of spatially 
heterogeneous dynamics (SHD) is evolving [6-14]. Sets of neighboring particles move with 
enhanced or diminished mobility relative to the average on a time scale intermediate between 
ballistic and diffusive motion. While there has been interest in the possible relation between 
clusters obtained from a SHD analysis and the CRR of the AG theory, a link between the 
quantitative SHD methods and the AG predictions has not been found. 

Here we show that on the time scale where SHD is prominent, the average cluster size 
n* can be related to the mass z of CRR, thus connecting the quantitative SHD analysis to 
the qualitative approach of AG. Our results are based on molecular dynamics simulations 
of the extended simple point charge (SPC/E) model of water [15]. We simulate N = 1728 
molecules at fixed density p = 1.0 g/cm 3 and a range of T from 200 K to 260 K (at 10 K 
intervals). For each T, we run two independent simulations to improve statistics. We find 
that D can be fit with 

D~(T-T MCT y (lb) 

using the values for the mode coupling temperature Tmct = 193 K and the diffusivity 
exponent 7 = 2.80 reported in Ref. [16]. 

To facilitate comparison with previous work, we use the same approach to define SHD 
clusters as that employed to study a Lennard- Jones (LJ) mixture [8] and experiments on 
colloids [14]. We define the mobility of a molecule at a given time t as the maximum 
displacement of the oxygen atom in the interval [to, to + At]. Following [8], we calculate the 
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self part of the time- dependent van Hove correlation function [17] G s (r, t) at t = t*, the time 
at which the non- Gaussian parameter 

Mt) = l(At))/(r\t)) 2 - 1 (2) 

has a maximum [18]. We fit G s (r,t) with a Gaussian approximation G (r) and define r* as 
the second intersection between these distributions. We find r* is in the range 0.2 — 0.25 nm 
for all T. We focus on the fraction of "mobile" molecules given by = 4irr 2 G s (r, t*)dr, 
i.e. the average fraction of molecules with a displacement larger than r* in the interval t* . 
Depending on T, we find 6% < < 8%. For simplicity we fix = 7% for all T. Similar 
values of were found in atomic systems [7,8,14] and in polymer melts [19]. Finally, we 
define a cluster of mobile molecules at each interval At as those mobile molecules whose 
nearest neighbor oxygen-oxygen distance at time to is l ess than the first minima of the 
oxygen-oxygen radial distribution function [20,21]. 

In Fig. 1 we present four snapshots of mobile particle clusters at T = 210 K for At = t*. 
particles follow each other in a roughly linearly fashion [7,14]. In LJ systems [7], mono atomic 
liquids [22] and polymers [23] complex clusters are composed of more elementary "strings" , 
where particles follow each other in a roughly linearly fashion. For water, small clusters are 
indeed string-like (e.g. Fig. 1(a)), but the molecules conform to the hydrogen bond geometry, 
and hence the clusters are generally less linear than clusters found in LJ systems. Like in the 
simpler systems, clusters become less string-like as their size increases [7,24], and the fraction 
of branching points — molecules with more than two neighbors (Fig. 1(b)) — increases with 
increasing cluster size. Indeed, the larger the cluster, the more complicated is its structure- 
becoming increasingly ramified (Fig. 1(c)) and even exhibiting loops (Fig. 1(d)) [25]. Hence 
it appears that the basic features of SHD found for models of simple liquids extend to the 
more complex molecular liquid, water. 

To relate SHD to the AG approach, we calculate the average cluster mass (ra(At)) for each 
T. In the AG approach to dynamics, the CRR are characterized by the number of particles 
z and configurational entropy s coll f(z) of the CRR; AG argue that z = Ns con f(z)/S con {. 
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Motivated by the recent finding that the average instantaneous cluster mass scales inversely 
with the entropy in a model of living polymerization [26], we use n* = (n(t*)) as a measure 
of z, since at t* correlations are very pronounced and (n(t)) is nearly maximal [27]. We find 
a linear relationship between n* and l/S con { (Fig. 2(a)), 

n* ~ 1 oc -L. (3) 

*Jconf 

This suggests that n* — 1 can be understood as a measure of z and provides a quantitative 
connection [28] between SHD clusters and the AG approach. It is necessary to subtract 
one from n* to obtain direct proportionality, implying that a cluster of unit size does not 
correspond to a CRR [8]. Equation (3) provides a clear link between a cluster property, n*, 
and a property of the PES, S con {. Furthermore, given that S con { and the diffusion constant 
D have been previously related [4], it follows that 

D ~ e - A ^-^ T . (4) 

Our results in Fig. 2(b) confirm this expectation. 

We next address the question of how SHD clusters in water depend on the observation 
time At. We focus on the number average (n(At)) and {n(At)) w = (n 2 (At))/(n(At)); 
(n(At)) w , the weight average cluster size, is the average size of a cluster to which a randomly 
chosen molecule belongs. Figure 3 shows (n(At)) and {n(At)) w for T = 210 K. To eliminate 
the random contribution, we normalize (n(At)) w by (n r ) w , the weight average cluster size 
for <f)N randomly chosen molecules. For comparison, we also include the non-Gaussian 
parameter a 2 (At) and the mean-squared displacement (r 2 (At)) (Fig. 3(a)) which displays 
the three characteristic time regimes, ballistic, cage and diffusive. 

The behavior of (n(At)) w / (n r ) w is analogous to that for polymer systems [19], with 
the exception that there is a clear increase in (n(At)) w /(n r ) w at the time scale on which 
molecules go from the ballistic to the cage regime. This additional feature is likely due to 
strong correlations in the vibrational motion of the first-neighbor molecules, owing to the 
presence of hydrogen bonds. In Fig. 3(c) we show (n(At)) w / (n r ) w for all T. For T < 240 K, 
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the maximum in (n(At)) w / (n r ) w increases in magnitude and shifts to larger time scales with 
decreasing T. The plateau at the crossover from the ballistic regime is nearly T-independent, 
as expected since the mean collision time is nearly T-independent. For T > 250 K, the 
maximum and the plateau merge, and hence it is not possible to separately distinguish 
these features. 

We find that t* is slightly larger than t max , the time where the maximum of (n(At)) w 
occurs. Both characteristic times correspond to the late-/3/early-a time regime of the mode 
coupling theory (MCT). We find (Fig. 4) 

t* ~ (T-Tmct)"' 5 , [5 = 2.7 ±0.1] (5a) 

and 

tmax- (T-Tmct)"*, [a; = 2.7 ±0.1] (5b) 

For the LJ polymer melt, simulations show that x = l/2a, where a is the scaling exponent 
predicted by MCT for the (3 time scale, suggesting that t max may be a measure of the 
"elusive" (3 relaxation time scale [19]. MCT predicts that knowing 7, defined in Eq. (lb), is 
sufficient for determining a [29]. If t max were a measure of the /3-time scale, then we would 
expect a = 0.28 [16], and hence x would be equal to 1.78. From Eq. (5b), we see that 
x > 1.78, so t max does not provide a measure of the (3 time scale. Additional tests of the 
temperature scaling of t max for other liquids are needed to determine the range of liquids for 
which t max can be considered a measure of the (3 time scale. 

In summary, the relation we find between n* and S con f, Eq. (3), provides a link between 
SHD and properties of the PES [30]. In the context of AG theory, our findings support 
the interesting possibility that n* is a measure of the size of the cooperatively rearranging 
regions. Furthermore, our simulations show that SHD in water are qualitatively similar to 
those found in LJ systems, but the cluster shapes are strongly influenced by the geometry 
of the hydrogen bond network. 
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FIGURES 




FIG. 1. Four clusters of mobile molecules found at T = 210 K defined with an observation 
time At = t*. Tubes connect neighboring molecules whose oxygen-oxygen distance is less than 
0.315 nm, the first minimum in the oxygen-oxygen radial distribution function. Small clusters can 
be either (a) string-like or (b) non-string-like, showing branching points (molecules with more than 
two neighbors). (c,d) Larger clusters exhibit more complicated structures. Clusters can usually be 
decomposed in substrings, although this can be complicated by the presence of loops. 
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FIG. 2. (a) The average cluster size n* is proportional to the inverse of the configurational 
entropy 5 con f suggesting that n* — 1 can be used as a measure of the size of the cooperatively 
rearranging regions hypothesized by Adam and Gibbs. (b) Log-linear plot of (n*— 1) /T as a function 
of the diffusion constant D. The AG prediction D ~ exp (A/TS con {) implies that logD ~ (n* — 1)/T. 
This relationship holds for almost three decades in D. 
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FIG. 3. (a) Mean square displacement (r 2 (At)) at T = 210 K showing the ballistic, cage 
and diffusive regimes (separated by dotted lines), (b) Average number of molecules (n(At)) (O) 
and normalized weight cluster size (n(At)) w (□). The behavior of all quantities correlate with 
(r 2 (At)). The maxima of (n(At)) w and (n(At)) occur at times slightly smaller than the time for 
the maximum in oc2(At) (CO) the non-Gaussian parameter, (c) Weight cluster size (n(At)) w / (n r ) w 
for temperatures ranging from 200 K to 260 K in intervals of 10 K. Note the T-independent plateau 
at the crossover from ballistic motion. 
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FIG. 4. Temperature dependence of t* and t maX ) the times at which the maxima of the 
non-Gaussian parameter and the weight average cluster size occur. We find t* ~ (T/Tmct — l) 5 
and i max ~ (T/Tmct — l) x where x = 5 = 2.7 ± 0.1. The expected value of the exponent x, 
following arguments in [19], is 1.78 (long dashed line). The values of t* and i max for T = 200 K 
were not included in the fits because 200 K is too close to Tmct = 193 K (deviations from MCT 
are known to occur as T — > Tmct)- 
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